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Abstract 

We study finite-volume effects on the masses of the ground-state octet baryons using covariant baryon 
chiral perturbation theory (ChPT) up to next-to-leading order by analyzing the latest nj = 2 + 1 lattice 
Quantum ChromoDynamics (LQCD) results from the NPLQCD collaboration. Contributions of virtual 
decuplet baryons are taken into account using the "consistent" coupling scheme. We compare our results 
with those obtained from heavy baryon ChPT and show that, although both approaches can describe well the 
lattice data, the underlying physics is different: In HBChPT, virtual decuplet baryons play a more important 
role than they do in covariant ChPT. This is because the virtual octet baryon contributions to finite-volume 
corrections are larger in covariant ChPT than in HBChPT, while the contributions of intermediate decuplet 
baryons are smaller, because of relativistic effects. We observe that for the octet baryon masses, at fixed 
m x L (S> 1) finite-volume corrections decrease as approaches its physical value, provided that the 
strange quark mass is at or close to its physical value, as in most LQCD setups. 

PACS numbers: 12.39.Fe Chiral Lagrangians, 12.38.Gc Lattice QCD calculations, 14. 20. Gk Baryon resonances 
(S=C=B=0), 14.20. Jn Hyperons 
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I. INTRODUCTION 



Understanding the origin of the masses of light hadrons has long been a central topic in strong- 
interaction physics. Due to the non-perturbative nature of the strong interaction at low energies, 
calculations based on first principles have only become possible with the advent of lattice Quantum 
ChromoDynamics (LQCD). LQCD has made remarkable progress in studies of strong-interaction 
physics in the past decade (see, e.g., Refs. tlL|2|]). Nowadays, fully dynamical calculations have be- 
come standard and therefore one of the most-difficult-to-estimate uncertainties related to "quench- 
ing" effects in LQCD calculations of early times have been removed. Nonetheless, LQCD sim- 
ulations still have to adopt unphysical simulation parameters: larger than physical light quark 
masses m u / d , finite lattice volume V = TL 3 , finite lattice spacing a, etc. To obtain physical re- 
sults, extrapolations to the physical point in terms of m u / d , L (T), and a must be performed!] i.e., 
m u /d -> m u/d (phys.), L(T) ->■ oo, and a 0. 

The extrapolation in light quark masses m u / d is usually termed as "chiral extrapolation" [see 
Refs. [|3|-l6l] for some early studies of the nucleon mass in SU(2)]. In the real world, chiral sym- 
metry and its breaking pattern govern the dynamics of low-energy strong interaction. This is 
systematically and consistently formulated in an effective field theory called chiral perturbation 
theory (ChPT) Ir7l-ll4ll. At present only a few LQCD simulations have been preformed directly at 



the physical light quark masses [|15L ll6ll. while most calculations still require some kind of chiral 
extrapolation, often introducing sizable uncertainties to the final results. 

LQCD simulations, by definition, are performed in a hypercube with its volume as a simulation 
parameter. The volume has to be large enough such that physics in a finite hypercube is approxi- 
mately the same as that in infinite space-time. To have effects of this origin under control, a rule 
of thumbj is that m n L should be larger than ~ 4. In the p-regime (where m n L ^> 1), it was 
first suggested by Gasser and Leutwyler that one could use ChPT to evaluate finite- volume correc- 



tions IU8L 



19|]. An alternative approach is the Liischer formula [|2C ] and its resummed version (for 



a recent reference in the context of the nucleon mass, see Ref. ||2l|l ). 

In recent years, it is found that three-flavor (u, d, and s) ChPT at next-to-leading order (NLO) 
has difficulties to accommodate recent LQCD results, particularly in the one-baryon sector. ^1 In 
the case of light hadron masses, it was shown that NLO heavy baryon (HB) ChPT cannot describe 

1 In this work, we limit our discussions to LQCD simulations of zero-temperature physics 

2 We will see that for the ground-state octet baryon masses studied in this work, this requirement depends on the value 



of m^: at fixed m^L the larger the m^, the larger the finite-volume corrections, as previously noted in Ref. 11711 

3 For an update on the present situation in the mesonic sector, see Ref. 02211 . 
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the latest LHP [23J and PACS-CS Q24J] lattice data. On the other hand, co variant baryon ChPT 
supplemented with the extended-on-mass- shell (EOMS) scheme Il25ll and ChPT regularized by 
a cutoff (long range regularization) [26] are shown to be able to describe much better the same 
lattice data. 

In terms of chiral extrapolation, it seems that the advantage of covariant ChPT over HBChPT 
in the one-baryon sector has been established, particularly in the three-flavor case. However, a 
detailed study of finite-volume effects using three-flavor baryon ChPT is still missing. To perform 
such a study, it is advantageous to have LQCD simulations performed with the same setup except 
for the lattice size (volume). Such LQCD results have recently been provided by the NPLQCD 
collaboration JTvIl , where simulations are performed with n f = 2 + 1 clover fermions in four lattice 
volumes, with spatial extent L ~ 2.0, 2.5, 3.0, and 3.9 fm, using an anisotropic lattice spacing of 
b s ~ 0.123 fm in the spatial direction and b t = b s /3.5 in the time direction, and at a pion mass of 
m 390 MeV. With these results, we can now perform a detailed study of finite-volume effects 
on the masses of ground-state octet baryons using both the covariant formulation of ChPT and its 
non-relativistic counterpart (HBChPT). 

In Ref. H27Q , a covariant formulation of ChPT using the infrared scheme [|28|l up to next-to-next- 
to-leading order (NNLO) was used to study the dependence of nucleon masses on the lattice size 
L in the n/ = 2 case. The authors concluded that NNLO relativistic ChPT can describe well the 
finite-volume effects of nucleon masses. In this work, however, contributions of the intermediate 
A(1232) were not considered. In Ref. Q29Q, the effects of the intermediate A(1232) were studied 
and the authors pointed out that its effects are important but can be encoded into the relevant low- 
energy constants and therefore explicit inclusion of the A(1232) into NNLO ChPT in the infrared 
scheme is not necessaryO 

In Ref. II 1711 . on the other hand, it was shown that three-flavor HBChPT at NLO can describe 
the observed volume dependences reasonably well and NLO HBChPT with the decuplet states 
integrated out does not provide a reliable description of the finite-volume effects at a pion mass 
of 390 MeV. We will show that indeed at NLO virtual decuplet baryons play a more important 
role in HBChPT than in covariant ChPT. In fact, in HBChPT much of the observed finite-volume 
effects can only be explained when contributions of virtual decuplet baryons are taken into account. 
This is because the virtual octet baryon contributions to finite-volume corrections are larger in 

4 The additional low-energy constant ej[ , introduced to compensate the logarithmic dependence on the renormaliza- 



tion scale, makes a direct comparison of their NLO results with those presented in this work difficult. 
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FIG. 1. Feynman diagrams contributing to finite- volume effects on the masses of the ground-state octet 
baryons up to next-to-leading-order. Solid lines denote octet baryons, solid double lines decuplet baryons, 
and dashed lines represent pseudoscalar mesons. 

covariant ChPT than in HBChPT, while the contributions of intermediate decuplet baryons are 
smaller, because of relativistic effects. 

This paper is organized as follows. In Section 2, we calculate the finite-volume effects on the 
masses of the ground-state octet baryons in the covariant formulation of baryon ChPT. In Section 
3, we compare finite-volume effects predicted by covariant baryon ChPT and those by HBChPT 
and study the latest NPLQCD results. We conclude in Section 4. 



II. FINITE-VOLUME CORRECTIONS TO GROUND-STATE OCTET BARYON MASSES 



In Ref. Il25h . the lowest-lying octet and decuplet baryon masses are calculated in a covariant 



formulation of ChPT supplemented with the Extended-on-Mass-Shell (EOMS) renormalization 

nn n 

scheme pQ . I31Q to conserve proper chiral power countingO It is shown that at NLO covariant 
ChPT can better describe the LQCD simulations compared to HBChPT^ We will not repeat the 



same calculation here and refer interested readers to Ref. Il25h for the details. In this work we 
concentrate on finite- volume corrections and spell out details of the calculation which have been 
skipped in Ref. Ir25ll . 

Physically, finite-volume corrections can be easily understood: Because of the existence 
of space-time boundaries, the allowed momenta of virtual particles become discretized, i.e, 

oo 

one has to replace a momentum integral of the form J dk by an infinite sum of discretized 

— oo 

oo 

momenta Yl (x) n ( assumm g periodical boundary conditions). In LQCD simulations of zero- 

n=— oo 

temperature physics the temporal extent is generally larger than the spacial extent such that the 

5 In this work both EOMS ChPT and covariant ChPT will be used to refer to ChPT supplemented with the EOMS 

prescription, unless otherwise explicitly specified. 

6 In the study of light quark mass dependence of the D and D s decay constants, it is also observed that covariant 

ChPT does a better job compared to the heavy-meson ChPT 13211 . 
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integral in the temporal dimension can be treated as if it extends from — oo and oo. As a result, 
only the integral in the spacial dimensions should be replaced by an infinite sum. Obviously only 
loop diagrams are affected by the existence of space-time boundaries. In studying the octet baryon 
masses up to NLO, only two such loop diagrams contribute and are shown in Fig.Q] 

To calculate the finite-volume corrections originating from these loop diagrams, one simply 
follows the conventional way of calculating Feynman diagrams and needs only to pay attention to 
the fact that now the temporal and the spacial dimensions must be treated differently. Diagrams 
(a) and (b) of Fig. Q] yield, generically, 

d A k ^-f + M^ 

(2tt) 4 (k 2 -m 2 M + ie)((p - k) 2 - M 2 + ie 



G N -i I 77m— 2 , TaTJZ u\2 n/ r2 . ■ "> C 1 ) 



r d A k j mae (p - k) m k e S ab (p - k)j nbf {p - k) n y 
D ~ l J (2tt) 4 (k 2 -m? M + ie)((p-k) 2 -M*+ie) () 
with 7 mae = I( 7 m 7 a 7 e - 7 e 7 a 7 m ), 7 " fe / = I( 7 n 7 V-7 / 7V),and S ab (p) = -(f + M D )(r] ah - 
_L_yy _ _L_ 2!faV _ £=|2g£) J23O. In Eqs. CEJ), M and M D are the octet and de- 
cupet baryon masses at the chiral limit, and m M is the mass of a Nambu-Goldstone boson. 



As in Ref. I125h . we have adopted the "consistent" coupling scheme advocated by Pascalutsa et 



al.M 



3511 to describe the interactions between octet and decuplet baryons. In infinite space-time, 



the above integrals have been calculated in Ref. |25|] and the results can be found there. As ex- 



plained in Ref. |25J], the above loop functions contain power-counting-breaking (PCB) terms and 
therefore additional steps need to be taken to conserve a proper chiral power-counting scheme. 
Among the different approaches, the EOMS scheme has been shown to be superior to the heavy- 



baryon or infrared approaches (see Refs. yl, for an in-depth discussion on this topic). 

In the following, G^^ 1 ^ in the covariant framework represent the loop functions in which 
PCB terms have been removed using the EOMS prescription. Instead of calculating the integrals, 
Eqs. (I1I2|) . in a finite hypercube, we calculate the following differences: 



5G N = G N (L)-G N (oo), (3) 

5G D = G D (L)-G D (oo), (4) 

where G N / D (L) and G N / D (oo) denote the integrals calculated in a finite hypercube and in infinite 
space-time. These quantities have several features that make calculations more feasible than a 
direct computation of Gn/d(L). First, because Gm/d{L) and Gjvm(oo) have the same ultravi- 
olet behavior, 5G N / D are finite and can therefore be calculated in four dimensions. Second, the 



unwelcome PCB terms appearing in a covariant baryon ChPT calculation are absent because they 
emerge from short-distance physics while such short-distance properties are the same in G n j d{^L) 
and G N / D (oo). As a result, PCB terms vanish in the differences 5G N / D and no power-counting- 
restoration schemes, such as EOMS or IR, are needed to calculate 5Gn/d- 
For G N , one has 



G 



N 



i I dx 



d 4 k 



f + M )lf; 



(2tt) 4 {(k - px) 2 - M 2 N ) 2 



i I dx 



d 4 k ${k 2 -2k-p) + 2k 2 M 
(2vr) 4 ((k - px) 2 - M 2 N ) 2 ' ( } 



where M. 2 N = x 2 Mq + (1 — x)m 2 M — it. Calculating the integral [Eq. ©] in a finite hypercube 
requires treating the temporal and spacial dimensions differently. We choose to work in the baryon 
rest frame, i.e., p^ = (M , 0). In this frame, 



G 



N 



i dx 



dk f dk (~f k -7-k)(k 2 -k 2 -2k M ) + 2(k 2 -k 2 )M 



2tt J (2tt) 



((k Q -xM ) 2 -k 2 -M 2 N ) 2 



(6) 



This can be easily calculated by performing a shift in k (k — >■ k' Q + xM ), Wick rotating k' Q 
(k' — > ik' ), and then performing the integration over k' . The result is: 



G 



N 



dx 



o 



dk 



1/2 



-m (2x + r 



k 2 + M 2 N 



M (M 2 x 6 + M 2 N (x + 2)) 



1 



k 2 + M 2 N 
(7) 



From this and utilizing the master formula as provided, e.g., in Ref. 113911 . one can easily obtain 
i 



6G 



N 



dx 



^M (2x + 1)S 1/2 (M%) - -\m (M 2 x 3 + M 2 N (x + 2))5 3/2 (M 2 N ) 



(8) 



where 



S r {M' 



2- 1 / 2 - r (VM 2 f- 2r 



71 



3 / 2 r(r) 



^(L\OW\n\)~ 3/2+r K 3 / 2 _ r (L\OW\n\), 



(9) 



oo oo oo 

where K n (z) is the modified Bessel function of the second kind, and = 2~2 2~2 S (1 — 

n^O n x =— oo n y =— oo n z =— oo 

S(\n\,0)) withn = (n x ,n y ,n z ). 



In a similar way, one can calculate 5Gr>: 
M 2 {M {1 -x) + M D 



5G 



D 



dx 



6M£ 



M 2 (M (l-x)+M D )M 2 D 



QM 2 D 



(10) 
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FIG. 2. Finite-volume corrections, <5G/v (left) and 5Gd (right), as functions of lattice size L for = 0.2, 
0.4 and 0.6 GeV. In the evaluation of 5Gd/n i n covariant ChPT Mq = 0.8 GeV. The decuplet-octet mass 
gap is set to be A = 0.231 GeV. 



where M 2 



D 



X 



l M 2 - x(M 2 - M 2 D ) + (1 - x)m 2 M 



it. 



The corresponding results in HBChPT can be found in Ref. [39] and are given below 



5G 



(HB) 
N/D 



c 



N/D 



d\/3 A J2 [(£|n|r^i(£/?A|n|) - /3 A K (L/3 A \n\)] 

n^0 



(ID 



M 



D 



- M Q for <5Gg IB) and A = for 5Gf B \ C N = -4, 



where /3 A = X 2 + 2AA + m 2 M with A 
andC D = -16/6. 

Numerically one can easily check that 5Gn/d [Eqs. (|8I10I) 1 turn out to be identical to SG^^ 
in the limit of M — > oo, which confirms our statement that there are no PCB terms in SG^/d- 

1 oo 

Furthermore, we have noted that changing the integration region in Eqs. (|8I10I) from J to J, 

o o 

which corresponds to the infrared prescription Il27ll . has negligible effects on the numerical results 
of those integrals. 



III. RESULTS AND DISCUSSION 



A. Finite-volume corrections: HB vs. EOMS ChPT 



Before looking at the NPLQCD results, it is instructive to compare the finite- volume corrections 
computed in covariant ChPT and HBChPT. In Fig.[2l finite-volume corrections 5Gn/d for different 
m w are plotted as functions of the lattice size L. The well-known rapid decrease of finite-volume 
corrections with increasing L is clearly seen. The finite-volume corrections can be parameterized 
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FIG. 3. Finite-volume effects, 5Gn/Gn (left) and 8Gd/Gd (right), as functions of e mnL /(m n L) for 
m,r = 0.2, 0.4 and 0.6 GeV. In the evaluation of 5G D/N (G D/N ) in EOMS ChPT M = 0.8 GeV. The 
decuplet-octet mass gap is set to be A = 0.231 GeV. 

as a linear function of e~ mnL /(m^L), as indicated by explicit calculations for large L with both 
ChPT and the Luscher methodO It is clear that for 5Gn at fixed ro, and L the covariant results 
are larger than the HB results. The differences become smaller as decreases. On the other 
hand, for SGd the covariant results are smaller than the HB results and the differences become 
larger as decreases. It should be stressed that the enhancement of virtual octet contributions 
and reduction of virtual decuplet contributions in EOMS compared to HB ChPT have important 
consequences for extraction of the decuplet-octet axial couplings from finite- volume dependence 
of octet baryon masses, as we will show in the following subsection. 

In ChPT, where finite- volume corrections are calculated order by order from loop diagrams, the 
ratios of 5Gn/d/Gn/d demonstrate better the differences between various ways of calculating the 
loop diagrams, e.g., HB vs. EOMS. In Fig.|3l we plot 5Gn/d/Gn/d as functions of e~ m7lL / (m^L) 
for = 0.2, 0.4, and 0.6 GeV, respectively. The finite-volume effects decrease almost linearly in 
both cases, reflecting the fact that to a large extent they can be parameterized as a linear function in 
terms of e ~ mnL / (m n L). Such a feature has been employed in Ref. Ill7ll to perform a phenomeno- 
logical fit to their data. There are, however, some subtle differences between the EOMS and HB 
results, which one cannot see easily from Fig. |2l At fixed m n L, in the covariant case the larger 
the m w , the larger the finite- volume effects. On the other hand, in HBChPT 5Gn/Gn depends 



7 To factor out such an exponential dependence, in the rest of this work, we will plot relevant quantities as functions 

of e- m ' L /{ m7r L). 
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FIG. 4. The NPLQCD octet mass data [17] fitted with NLO covariant ChPT and HBChPT. The bands are 
the full results at the 68% confidence level and the solid (dashed) lines are the best fits with C = 0. 

only on m^L^] while the dependences of 5Gd/Gd on rn n are clearly different from its covariant 
counterpart. 



B. Study of the NPLQCD results 

Now let us turn to the NPLQCD data. The lattice results in units of temporary lattice spacing 
are given in Table I of Ref. II 1711 . which can be translated into physical units using b t = b s /3.5 
with b s = 0.1227 ± 0.0008 fm. We have chosen to fit the masses with physical units, which is 
equivalent to fitting the masses with lattice units. Of course, to get the final physical results one 
must take into account also the uncertainties due to the determination of the lattice spacing. This, 
however, does not affect the quality of our fits. 

The octet baryon masses at NLO covariant baryon ChPT in infinite space-time have the follow- 
8 This can be immediately understood by looking at the corresponding analytical results: 

G^ B) = -2nml, ( 12 ) 
«5G^ B) = 2nml ^{Lm^y 1 exp(-m 7r i|n|). (13) 
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ing form: 

=mo- £ C^l + T^ £ ^^K), (14) 

a=b,c 

where £g a ^ and are tabulated in Table V of Ref. [25], the loop functions are given in 
Eqs. (A3, A4) in the same reference, which contain PCB terms that have to be removed by the 
EOMS prescription as explained in detail in Ref. 12511 . To calculate the baryon masses in a finite 
hypercube, one simply replaces the loop functions H^' with their counterparts calculated in a 
finite hypercube as provided in the previous section, i.e., 

H®^B® + \8G N {L), (15) 

H®->H® + 1SG D {L). (16) 

To fit the octet baryon masses in the covariant ChPT one has four low-energy constants (LECs) 
to determine: bo, bo, bp, M . The other parameters are given the same values as those used in 
Ref. [25] , i.e., D = 0.8, F = 0.46, fi = 1 GeV, F w = F K = F v = 1.17/ w with f v = 92.4 MeV@ 
Since in the present case all the lattice data are obtained at the same pion mass, one could not 
distinguish bo and M and therefore we have taken M = 0.8 GeV, as suggested by the covariant 



ChPT study of the PACS and LHP data 1251140(1. A moderate variation of M from its central value 
by, e.g., 0.2 GeV, does not change the results in any appreciable way. It should be noted that M 
does not enter loop calculations in HBChPT. For the decuplet-octet mass gap, we take the average 
value A = 0.231 GeV. Using the mass gap between the nucleon and the A(1232), 0.291GeV, for A 
has minor effects on our fits. Furthermore, to study the effects of the intermediate decuplet baryons 
we allow C, the decuplet-octet axial coupling, to vary. At the end, we have four parameters to fit 
16 lattice data. 

The fitted results using both covariant ChPT and HBChPT are shown in Fig. HI Both methods 
provide a reasonable fit to the lattice data with similar quality, compatible with the phenomeno- 
logical fit performed in Ref. [|17l]. though with less parameters (four vs. eight). On the other hand, 
for the E, ChPT seems to predict a larger finite-volume dependence than suggested by the lattice 
data. The values of the LECs corresponding to our best fits are given in Table I. Two things are 
noteworthy. First, the values of b , b D , and bp are close to those obtained in Ref. Ii25ll . Second, 



9 It should be noted that in the present study, SU(3) breaking effects are entirely induced by the masses of the pseu- 
doscalar mesons, while for the couplings we have not introduced any explicit SU(3) symmetry breaking, because 
they can not be determined by the NPLQCD data. 
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TABLE I. Values of the low-energy constants from the best fit to the NPLQD data with x 2 /d.o.f . « 1.6. 





bo 


bo 


b F 


C 


EOMS 


-0.81 ±0.08 


0.12 ±0.04 


-0.46 ± 0.04 


0.73 ±0.23 


HB 


-1.85 ±0.20 


0.61 ±0.11 


-1.01 ±0.09 


0.94 ±0.09 



the central value of C from the best fit in EOMS ChPT is only about 3/4 of the value we used in 
Ref. [25], which is fixed from the A decay width. Because C should be understood as an average 
of all the decuplet-octet axial couplings, such a value, taking into account its uncertainty, is not 
out of the range of our expectations. On the other hand, the value of C fixed from the HPChPT 
fit is slightly larger, 0.94 ± 0.09. This is larger than the value used in most HB calculations, e.g., 
C = gAN/% = 0.7 [39] and C = 0.76 [26]. Both results are, however, roughly consistent with 
each other, keeping in mind that NNLO ChPT contributions could be as large as 30% of those of 
the NLO. 

The effects of the virtual decuplet baryons can be best seen by fitting the NPLQCD data with 
C = 0. The corresponding results are shown by the solid (dashed) lines in Fig. @] It is clear that 
in the fit the octet-decuplet transition plays a larger role in HB than in covariant ChPT. In fact, 
in HBChPT virtual decuplet baryons play an even larger role than those of virtual octet baryons, 
which seems to be a bit unnatural (for a relevant discussion, see, e.g., Ref. Bill ). In Ref. II 1711 . it 
was concluded that the decuplet contributions must be taken into account. Our studies show that 



TABLE II. Extrapolated octet baryon masses in the infinite lattice size limit. 





EOMS ChPT 


HBChPT 


SU2 HBChPT [17] 


M N 


1152.7 ±1.4 


1151.7 ±1.5 


1151.3 ± 1.1 


M A 


1242.7 ±1.0 


1244.5 ± 1.0 


1241.9 ±0.8 


M s 


1278.2 ± 1.4 


1279.6 ± 1.4 


1280.3 ± 1.0 




1349.8 ± 1.2 


1349.5 ± 1.2 


1349.6 ± 0.7 
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FIG. 5. Predicted finite-volume corrections for the ground-state octet baryons using the best fitted LECs 
determined from the NPLQCD data 111 711 . where m s is assumed to have its physical value, and itlk (m v ) 
are related to through leading-order ChPT: m 2 K = \ (2Bom s + m 2 ), m 2 = ^{Am 2 K — m 2 ), and 
2B m s = (2m 2 K - m^)| phys . 

this is indeed the case, but more so in the HBChPT than in the covariant ChPT. From the above 
discussion, we reach the same conclusion as Ref. [170 that extraction of the decuplet-octet axial 
couplings from the present lattice data by studying the volume dependence of the octet baryons 
masses cannot be taken too seriously. 

The extrapolated octet baryon masses in the limit of L — > oo are shown in Table II. For the 



sake of comparison, we have tabulated the results of the NLO SU(2) HBChPT fit from Ref. QJ/ZD - 
The uncertainties are from the lattice data (systematic plus statistical) and we have ignored our 
systematic uncertainties and the uncertainties due to the lattice scale in converting the lattice results 
from lattice units to physical units[j°| We note that the extrapolated masses using different methods 
are roughly consistent with each other. This is mainly because the lattice simulations have been 



As can be seen from Table III of Ref. |l7D, such uncertainties are far more important numerically. 
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performed with a large volume, L = 3.9 fm, where finite-volume corrections are almost zero and 
which strongly constrain the extrapolations. 

In Fig. |51 we show the predicted finite-volume corrections for different using the LECs 
given in Table I. The HB and covariant ChPT predictions are quite similar, except for the £ and 
S, where the difference could reach 10 — 20% at a pion mass of 0.4 GeV and small L. Therefore, 
it can be concluded that at NLO one can use either HBChPT or covariant ChPT to describe finite- 
volume corrections to the octet baryon masses, keeping in mind that relativistic effects are large 
and one must be careful about interpretation of the extracted physical quantities. 

IV. SUMMARY AND CONCLUSIONS 

We have studied finite- volume effects on the octet baryon masses by analyzing the latest n/ = 
2 + 1 NPLQCD data with a covariant formulation of baryon chiral perturbation theory and with 
heavy baryon chiral perturbation theory. It was shown that although both approaches can describe 
the lattice data reasonably well, the underlying physics is different: Decuplet contributions play a 
less important role in covariant ChPT than in HBChPT at next-to-leading order because relativistic 
corrections enhance virtual octet contributions and reduce intermediate decuplet contributions. 
This makes it difficult to reliably extract the values of the decuplet-octet axial couplings from the 
volume dependences of the octet baryon masses. Simulations with a larger statistics, multiple 
volumes, and different pion masses will likely better serve such purposes. 

We have shown that at fixed m n L (3> 1) finite- volume corrections become smaller for smaller 
m w , as the NPLQCD collaboration has pointed out. For the LQCD calculation of the nucleon 
mass, m^L > 4 is needed to have an finite- volume correction at the order of 1% at a pion mass of 
~ 0.3 GeV, while at a pion mass of ~ 0.2 GeV, m^L > 3 is enough. For the other octet baryons, 
the dependence of finite- volume corrections on e ~ m7rL / (m^L) is weaker than that of the nucleon, 
provided that m s is close to its physical value. 

It should be stressed that the fact that both HBChPT and covariant ChPT can describe finite- 
volume effects on the ground-state octet baryon masses does not mean that both approaches are 
capable of describing the pion-mass dependence of these quantities. In fact, it has been shown that 
at NLO covariant ChPT is more suitable for the purpose of chiral extrapolation. 

At present, to go to higher chiral orders in the applications of three-flavor covariant baryon 
ChPT, one faces the problem of a large number of poorly known low-energy constants, such as 

13 



the decuplet-octet axial couplings. LQCD simulations with different quark masses and a sequence 
of different volumes help setting constraints on the values of these LECs as shown in the present 
work. 
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